Electron-Electron scattering and resistivity of ballistic multimode channels 
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We show that electron-electron scattering gives a positive contribution to the resistivity of bal- 
listic multimode wires whose width is much smaller than their length. This contribution is not 
exponentially small at low temperatures and therefore may be experimentally observable. It scales 
with temperature as for three-dimensional channels and as T^/^ for two-dimensional ones. 
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I. INTRODUCTION 



Electron-electron scattering does not contribute to 
the resistivity of a homogeneous conductor with a pa- 
rabolic spectrum. The reason is that the current den- 
sity is proportional to the total momentum of elec- 
tron gas, which is conserved in electron-electron colli- 
sions. If umklapp processes^ are neglected, this con- 
tribution may be nonzero only in a presence of mi- 
croscopic or macroscopic inhomogeneities. In particu- 
lar, quantum-mechanical interference between electron- 
electron and electron-impurity scattering results in a 
temperature-dependent correction to the conductivityi^i^ 
The electron-electron scattering is also known to affect 
the resistance of narrow channels with boundary scatter- 
ing because it deflects electrons moving along the chan- 
nel axis to the boundaries, where they can dissipate 
their momentum^ In short and wide ballistic contacts, 
the scattering in the electrodes also influences the resis- 
tance because the collisions change trajectories of elec- 
trons near an aperture in a diaphragm i^"— The electron- 
electron scattering is also known to influence the conduc- 
tance of quantum contacts with saddle-point potential 
sharp enough to violate momentum conservation)^ 

In all the above cases, the inhomogeneities are intro- 
duced into the system explicitly. However in finite-length 
conducting channels, the inhomogeneity is due to a mere 
presence of electron reservoirs at the ends. In^, Lunde et 
al. found that the contribution of electron-electron scat- 
tering to the conductance of a two-channel quantum wire 
attached to reservoirs becomes nonzero at some specific 
relations of the Fermi wave vectors of electrons in dif- 
ferent channels. More recently, a number of authorsiSrJ^ 
addressed the effect of electronic collisions on the conduc- 
tance of a single-channel quantum wire of a finite length. 
As two-electron collisions do not affect the current in a 
strictly one-dimensional system because of momentum 
and energy conservation, this contribution is determined 
by triple electronic collisions, which involve an empty 
electronic state near the bottom of the band. Therefore 
it is exponentially small at low temperatures and more- 
over, it is zero for some scattering potentials, e.g. for the 
point-like one. Importantly, this effect does not require 
an explicit presence of an inhomogeneity that can ab- 
sorb electronic momentum. Instead, this inhomogeneity 
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FIG. 1: A sketch of the ballistic multimode wire with two 
electrons colliding half-way. The collision of the left-moving 
electron (blue) with the right-moving one (red) prevents the 
left-mover from crossing the section of the wire shown by the 
dashed line in the middle and therefore affects the current 
flowing through it despite the momentum conservation. 



is implicitly introduced through the boundary conditions 
at the ends of the wire, which state that electrons coming 
to the wire from a reservoir have the same equilibrium 
distribution as in the reservoir. 

In this paper, we calculate the resistivity of a finite- 
length narrow multichannel ballistic conductor that re- 
sults from weak electron-electron scattering. Because the 
motion of electrons is also possible in the transverse di- 
rection, two-particle collisions affect the current despite 
the conservation laws. These collisions involve only elec- 
trons near the Fermi level, and therefore the effect is not 
exponentially small at low temperatures. 

The paper is organized as follows. In Sec. II we present 
the model and basic equations, in Sec. Ill we discuss the 
results, and Sec. IV presents the summary. Appendices 
A and B contain details of calculations for the 3D and 
2D cases. 



II. MODEL AND BASIC EQUATIONS 

Consider a metallic wire of uniform cross-section con- 
necting two massive electronic reservoirs. The length L 
of the wire is much larger than its transverse dimension 
a but much smaller than both the elastic and inelastic 
mean free paths. In its turn, a is much larger than the 
Fermi wavelength Xp. It is assumed that electrons are 
specularly reflected at the boundaries of the wire, so they 
cannot transfer their longitudinal momentum to the lat- 
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tice. In a two-dimensional (2D) case, such a wire could 
be realized by twisting up a 2D conducting strip to form 
a tube so that the boundary scattering could be elim- 
inated at all. We neglect here the effects of electron- 
electron scattering near the contacts of the wire with the 
reservoirs^. The kinetic equation for the distribution 
function of electrons inside the wire is of the forro^ 
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where the collision integral is given by 

Iee{v)=aeeVd ' ' ' 
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Qfee is the dimcnsionless interaction parameter, d = 2 
or 3 is the dimensionality of the system; 1/3 = mpi?/7r^ 
and V2 = Tnji^ arc the three- and two-dimensional two- 
spin electronic densities of states (ft = 1). The term 
with electric field is omitted because we address the linear 
response and the voltage drop is included in the boundary 
conditions 

/(p, >0,a; = 0) = /;^(ep-Sf^), (3) 

/(p, <0,X = i) = /i(£p-i?;^), (4) 

where x is the longitudinal coordinate, fni^s) = fo{s — 
eV) and /^(e) = /o(e) = 1/[1 + exp(e/r)]. The current 
through an arbitrary section of the conductor is given by 
an integral over the transverse coordinates 
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Equation ([T]) may be solved iteratively in the collision 
integral. In the zero approximation, the electrons move 
inside the wire along broken lines while retaining their 
energy and longitudinal momentum (see Fig. 1), so their 
distribution depends solely on the reservoir from which 
the trajectory of an electron with a given momentum 
originates. Therefore the coordinate-independent distri- 
bution is of the form 

/(°)(p,r) = e(p,)/fl(£p-SF) 

+ 0{-p,)fL{ep~EF) (6) 

and the zero-approximation two-spin conductance is just 
given by the Sliarvin forniul a^^'^^ 
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To the first approximation in aee, the correction to 
the distribution function 5f{p,r) at a given point r of 



the conductor is readily obtained by integrating Jee{/^°''} 
along the trajectory of an electron with momentum p 
that comes to point r from the corresponding reservoir. 
As Iee{f^^^} is coordinate- independent, this integration 
is reduced to multiplying it by the time of motion inside 
the conductor 



Sf{p.^ > 0,P±,X) ^—ieeipx > 0, (8) 

L — X 

(p, <0,K|). (9) 



Upon a substitution of Eqs. ([8]) and ^ into Eq. ^ 
Vx nicely drops off and taking into account the particle- 
number conservation by the collision integral, one obtains 
an expression similar to that ofi^ 
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which tells us that 61 is proportional to the rate of change 
in the number of electrons with px > 0. One might think 
that it contradicts the momentum-conservation law, but 
this is not the case. Though collisions between elec- 
trons do not change their total momentum, they change 
their trajectories and hence may prevent some of them 
from passing through the cross-section of the conductor 
at which the current is calculated (see Fig. 1). 

Equation ([5]) must be substituted into pUj) . Because 
of a piecewise distribution (jG)), it is convenient to present 
the result as a sum of integrals over left-moving and right- 
moving electronic states. Upon a decomposition of mo- 
mentum integration into the integrations over the angles 
and over the energies, one obtains 
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where the indices a, /3, and 7 correspond to left-moving 
(L) and right-moving (R) states. The distribution- 
dependent factor is given by 

FRap^{e, e'|ei, £2) 

= [1 -/«(£)] [1 -/„(£')] /Mei) A (£2) 

- fnie) fa{e') [1 - Me,)] [1 - ^(£2)] (12) 

with /fl(£) = fo{£-eV) and /l(£) = /o(£)- The quantity 

f d^p f d^k f d^p' 
Bn.^,{e,e |£i,£2) / / 

X J d^k' 5(£p -Ep-e) J(£k -Ef- e') 
X 6{e'^ ^Ef~ £1) 5{e'^ - Ef - £2)6(17,) e((7, fc,) 

X e(CT^p;)e(CT7fc;)(5(p + k-p'-k') (13) 
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with (Jn = 1 and ctl = —1 is the factor representing 
the phase space available for scattering. The quantity 
Brrll is zero because of momentum conservation, and 
the distribution-dependent factors Frrrji, Frlrl, and 
Frllr vanish because of energy conservation. Hence 
only four terms in the sum containing odd numbers 
of left-moving and right-moving states are nonzero. Note 
that Bajj-^s is symmetric with respect to permutations of 
the first and second pairs of arguments and to permuta- 
tions inside these pairs and also symmetric with respect 
to a simulataneous reversal of all its indices. Using these 
symmetry properties and linearizing FRapj with respect 
to the voltage, one may reduce the sum pT|) to a single 
term 
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where 



F{e,e'\ei,e2) = [l-fo{e)] [1 
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The quantity B = Brrrl exhibits different behaviors 
for d = 3 and d = 2. In the three-dimensional case, it is 
dominated by collisions of electrons with total momen- 
tum Q of the order oi pp and remains finite even if all 
the initial and final states lie exactly at the Fermi sur- 
face. This may be understood as follows. Consider two 
electrons with momenta p and k lying exactly on the 
Fermi sphere that scatter into p' and k', respectively. 
Because of the conservation laws, the two initial and the 
two final states form a rectangle inscribed in the sphere 
and lying in a plane perpendicular to Q, which is located 
at a distance Q/2 from its center (see Fig. 2). Depend- 
ing on the direction and length of Q, the plane Px ~ 
may cut the rectangle into two parts and isolate two, one, 
or none of the vertices from it. The quantity Brrrl is 
determined by momentum configurations with just one 
isolated vertex, which exist even in the zero-energy limit. 
To the lowest order in e, calculations give (see Appendix 
A) B = C3 m^pp where C3 = 0.222, and the relative 
correction to the conductance is 
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where Ipe is the standard length of electron-electron scat- 
tering. This means that a significant proportion of 
electron-electron collision contributes to the current. 

In the two-dimensional case, B vanishes at e = e' = 
ei = £2 = 0, and therefore the correction to the current 
is determined by collisions with small total momenta. 
Explicit calculations give (see Appendix B) 
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FIG. 2: A 3D Fermi surface and the four states involved in 
a collision of two electrons. States p, p' and k' correspond 
to electrons moving from left to right whereas state k corre- 
sponds to an opositely moving electron, (a) A 3D view of the 
Fermi sphere and (b) its section with the plane containing the 
four states with the trace of the plane px = (yellow line). 
The collision increases the number of right-movers by one. 



where 



0, 



+ y+^/v- \x\, y>\x\ 
\y\ < \x\ 



y, 



(18) 
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The difference from the three-dimensional case is readily 
seen in Fig. 3. All the four momenta may lie exactly at 
the Fermi circle only at Q = 0. In that case, the center 
of the rectangle coincides with the center of the Fermi 
circle and the line = cuts it into two equal parts 
so that two vertices lie in the half-plane p^ > and the 
other two, in the half-plane px < 0. If the four momenta 
are located within Se ^ T from the Fermi surface, a 
rectangle with one isolated vertex may be placed in this 
ring-shaped region only if all its vertices are located near 
the line Px = 0, where the portions of the ring at the 
opposite sides of the Fermi surface are almost parallel. 
Therefore characteristic values of Q that dominate B are 
determined by the length at which the bending of these 
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FIG. 3: A 2D Fermi surface and a scattering process that 
changes the number of right-moving electrons. The dashed 
areas show the parts of the surface that contribute to such 
processes. 



portions is of the order of their width Se/vp, i.e. Q ^ 
y/pps/vF- A substitution of Eqs. ([TT]) and dH]) into (fT4|) 
gives 

SG2 _ C2 L a,,T^/^ 

where C2 = 2.23. This suggests that at low tem- 
peratures, SG2/G02 is smaller than L/l^e by a factor 
(T/EfY^^ and only a small fraction of electron-electron 
collisions affects the current. 

The proportionality of SG3 and SG2 to L is the result 
of the linear approximation in a^e- In fairly long wires, 
this approximation breaks down. It is reasonable to be- 
lieve that at i ^> 00 the corrections to the conductance 
saturate to a finite value as is in the case of a single-mode 
channelii and therefore present a boundary effect. 



III. DISCUSSION 

In contrast to the positive contribution to the con- 
ductance from the scattering in short and wide ballistic 
contacts^"—, the negative contribution in long and nar- 
row wires docs not require an explicit transfer of electron 
momentum to a macroscopic obstacle. Instead, the ex- 
cess momentum is absorbed in the reservoirs just like 
the energy of injected electrons. The principal differ- 
ence from the case of short contacts is that the collisions 
in a long and narrow channel involve only two type of 
electrons, left-movers and riglit-movers,which both con- 
tribute to the current with opposite signs. In contrast 
to this, in wide and short contacts the collisions involve 



also a third type of electrons, by-passers. They by far 
outnumber the first two types, do not contribute to the 
current in the absence of scattering but can be converted 
into a left-mover or right-mover upon a collision. 

The scattering mechanism in wires is more efficient in 
3D than in 2D because the fraction of collisions contribut- 
ing to the current increases with the dimensionality. Un- 
like this, the scattering contribution in wide contacts is 
smaller in 3D than in 2D because the concentration of 
injected nonequilibrium electrons falls off faster in 3D 
reservoirs. 

The scattering contribution to the resistance may be 
well observable in GaAs heterostructures under realis- 
tic experimental condition. For an electronic density 
Us = 10^^ cm~^ and T = 4 K high-mobility samples 
may have the elastic mean free path as large as 20 fim. 
For a conducting channel of that length and aee ^ 1, 
the correction to the conductance should be about 10%. 
This correction can be isolated from the noninteracting 
Sharvin resistance ([7]) by changing the length of the chan- 
nel L in a 2D electron gas by means of electrostatic gates 
and isolating the component of the resistance propor- 
tional to L. 



IV. SUMMARY 

In summary, we have shown that two-electron colli- 
sions contribute to the electric resistivity of multimode 
finite-length metallic wires even in the absence of umk- 
lapp processes. The finite contribution to the resistiv- 
ity is due to the breaking of translational symmetry of 
the system by the presence of reservoirs. This effect has 
a power-law temperature dependence and is more pro- 
nounced in 3D wires than in 2D ones. 
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Appendix A: Calculation of B for the 3D case 

The phase space available for scattering may be 
rewritten in a form 

f d'^Q 

B{e,e'\£i,e2) = J j^;^ ARR{e,e' M) 

X AiiL(ei,£2,Q), (Al) 
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where Q = p + k = p' + k' is the total momentum of the 
coUiding electrons and 

A^p{s,e',Q) = J d'^p J d''keia^p,)eiap k^) 

X 5{ep -Ep-e) 5{e^, - Ep ~ e') 6{p + k~ Q). (A2) 
In 3D, explicit caleulations give 

2-Ve(A)/.,efg^4^i;±£-eos,l 



(27r)3Q 



B cos 05 — —— 



9 ? 
Pe - Pc 



X 9 cos ip - ^ , 



(A3) 



and 



^..(e,e',Q) = ^e(A) 

2 9 

X Pe - pi' 



X dipQ 



cos (p 



(A4) 



where = ^2m(e + -E_f), Pe' = y/2mJe^ + Ep), and 

A(g,£,e') = fQ'-bs-Pe')'" 

(Pe+P,0'-Q']- (A5) 

The calculations are greatly simplified if e — e' — )• 0. 
By going to spherical coordinates in (|A1[) and intro- 
ducing new integration variables u = Q/2pp and w ~ 
q2 /^-i/2 (,q^ g instead of the radial coordinate Q and the 
polar angle 6, one arrives at an expression B = C3 m^pF, 
where 
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, , arcsmi) arccosu „ „„„ / . „n 
X / dv — „ = 0.222. (A6) 



the form 



^i?ii(£,£ ,Q) = TTT' 
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27r2 ^A(g,e,e') 
X {e[(g2 + p2 _ p2^) + IQ^I ^A(Q,£,£') 
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(Q' -p' + IQyl VA(g, £,£')]}, (Bi) 
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^A(g,£i,e2) 

{e[(g2 g. + ig^i VA(g,ei,£2) 



X e 

+ 6 



X e 



-|gj VA(g,£i,£2)]}. (B2) 



The integral (|Aip is conveniently calculated by going to 
radial coordinates. As Eqs. ()Bip and (jB2[) depend on 
the \Qy\IQx ratio only through 0-functions, it takes up 
a form 



B(£,£'|£i,£2) 



2pF 



dQ- 



QS^iQ) 



= , (B3) 



16^' J '"'^ VA(Q,£,£')A(g,£i,£2 



where Sip{Q) is the angular width of the domain of in- 
tegration defined by the theta-functions in Eqs. (jBip 
and (jB2[) . The integral is dominated by small values of 



g ^ \/pf^I'Vf <C PFi so that 



Appendix B: Calculation of B for the 2D case » 1 i\ a i \ a 2 r\2 it, a\ 

A{Q, £, £') w A(g, £1, £2) w 4p|,g^ (B4) 

Equations (|Aip and (|A2|) are also valid in the 2D case, 
but the explicit expressions for Arr and Arl are now of The integration results in Eqs. (fT7|) and (|18p . 
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